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Aims. It is necessary to understand the dynamics of the atomic gas to use complex modeling and to carry out detailed compar- 
isons between theoretical models and observations. 
• Methods. In a companion paper, we present high resolution bidimensional numerical simulations of the interstellar atomic hy- 
O-i drogen. Here, we further characterize these simulations and we compare our results with various observations. 
Q Results. We give statistics of the column density and velocity along the line of sight and show that they compare favorably with 
5_l ■ observations of high-latitude lines of sight. We compute synthetic HI spectra and qualitatively discuss the information that could 
| be inferred if these spectra were observed. Finally, we extract CNM clouds and study their physical properties finding strong 
d . similarities with real clouds. In particular, we find that the clouds follow Larson-type relations, i.e M oc L 7 , where 7 ~ 1.7 (we 
^ speculate that in 3D, 7 ~ 2.5) and V < Sv 2 > oc L 0A . We also find that the distribution, Af(N), of the column density, N, of 
. , the CNM structures formed in the simulation follows Af(N) oc 7V~ 1,2 which is marginally compatible with the observational 
■ result obtained by Heiles & Troland (2005). From the mass-size relation and the mass spectrum, we derive an exponent for the 
?H ' column density distribution close to the value obtained in the numerical simulation. 

Conclusions. We conclude that the simulations reproduce various observational features reasonably well. An important implica- 
tion suggested by our results is that the "turbulence" within the cold interstellar atomic gas is mainly the result of individual 
long living cloudlet (confined by an external warm medium) motions rather than supersonic turbulence within nearly isothermal 
clouds. Another important aspect is that the CNM structures produced in the simulation present various physical characteristics 
that are similar to the characteristics of the molecular clouds. This raises the question as to whether the physical properties of 
the molecular clouds are determined at a very early stage, before the gas becomes molecular. 

Key words. Hydrodynamics - Instabilities - Interstellar medium: kinematics and dynamics - structure - clouds 



1. Introduction 

Since its first use in the 1950s (Ewen & Purcell, 1951) 
the 21 cm transition of atomic hydrogen (HI) has been 
extensively used (e.g. Crovisier 1981, Kulkarni & Heiles 
1987, Dickey & Lockman 1990, Joncas et al. 1992) to study 
the neutral atomic interstellar gas by numerous authors. 
Recent data have considerably improved our knowledge 
of this medium. Heiles & Troland (2003, 2005), using the 
Arecibo telescope, have done an extensive survey of HI 
clouds which has given reliable statistical results. To di- 
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rectly probe smaller scales of HI emission, various studies 
have been done using interferometers. Miville-Deschenes 
et al. (2003) and Taylor et al. (2003) use the Dominion 
Radio Astrophysical Observatory (DRAO) interferome- 
ter, Stil et al. (2006) use the Very Large Array (VLA) 
whereas McClure-Griffiths et al. (2005) use the Australian 
Telescope Compact Array (ATCA). 

The observations of HI have revealed its complex 
multi-phase nature. In HI, cold and dense (CNM) struc- 
tures are embedded in a warm (WNM) phase, with which 
there are approximately in pressure equilibrium, as pre- 
dicted by detailed computations of thermal balance (Field 
et al. 1969, Wolfire et al. 1995, 2003). The observa- 
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tions have also revealed that this multi-phase medium is 
strongly turbulent making their interpretations even more 
challenging and indispensable the use of adapted numeri- 
cal simulations. 

Numerous numerical studies of the large scale ISM 
(e.g. Vazquez-Semadeni et al. 1995, Rosen & Bregman 
1995, de Avillez & Breitschwerdt 2005a) and of the molec- 
ular clouds (e.g. MacLow & Klessen 2004) have been per- 
formed over the years and comparisons between observa- 
tions and simulations have been carried out both for large 
scale simulations (e.g. Ballesteros-Paredes et al. 2002, de 
Avillez & Breitschwerdt 2005b) and for molecular clouds 
(Padoan et al. 2003, Ossenkopf & Mac Low 2002). 

It is, however, only recently that numerical simulations 
devoted to the description of the warm and the dense 
atomic gas with a numerical resolution sufficient to de- 
scribe the physical scales of the problem (i.e. the CNM 
structures and the Field length) have been performed 
(Hennebelle & Perault 1999, 2000, Koyama & Inutsuka 
2000, 2002, Audit & Hennebelle 2005, hereafter paper I, 
Heitsch et al. 2005, 2006, Vazquez-Semadeni et al. 2006). 
For this reason, no tight comparison between observations 
and theory have yet been carried out for HI at scales below 

10 pc. 

In a companion paper (Hennebelle & Audit 2006, here 
after paper II), we present 2D high resolution numerical 
simulations aiming to describe self-consistently a turbu- 
lent atomic hydrogen flow from scales of few tens of parsec 
down to hundredth of parsec. The high resolution which is 
necessary to describe properly HI permits to obtain a good 
statistical description of the flow and to reach small scales 
close to the size of various small scale structures that have 
been observed in HI, like the so-called tiny small atomic 
structure (TSAS) and low column density clouds recently 
observed by Braun & Kanekar (2005) and Stanimirovic 
& Heiles (2005). Therefore, in spite of the fact that these 
simulations are only bidimensional, they constitute a good 
starting point to carry out preliminary comparisons with 
the available HI data. Note that the dynamics of spatial 
scales which is covered by these simulations is about ten 
times larger than what can be done in 3D yet. 

In this paper, we attempt to further characterize the 
results of the numerical simulations presented in paper 

11 and to perform comparisons between these simula- 
tions and the HI observations obtained in the Millennium 
Arecibo 21 centimeter absorption- line survey (Heiles & 
Troland 2003, 2005) since various interesting statistical in- 
formations have been extracted from these data. We also 
make a special effort to provide enough material to allow 
the observers to make easier the comparisons with their 
observations. In Sect. 2, we present some statistics of col- 
umn density and velocity along the lines of sight since 
this information is the most direct observable. In Sect. 3, 
we present various examples of synthetic atomic hydro- 
gen 21cm line spectra in emission and absorption, showing 
in each case the corresponding density and velocity fields 
along the line of sight. We discuss the influence of the 
convolution by a Gaussian beam of various width on the 
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Fig. 1. Top panel: column density, TV, along the lines of 
sight as a function of position. Bottom panel: integral of 
n/T along the line of sight. Full line represents lines of 
sight parallel to the y-axis, whereas dotted line represents 
lines of sight parallel to x-axis. 

synthetic spectra. Finally, we focus in Sect. 4 on the CNM 
structures formed in the simulations and their properties 
trying to make the link with some of the cloud properties 
available from the Millennium survey. Section 5 concludes 
the paper. 

2. Lines of sight 

In paper II, the structure of the flow which results from the 
numerical simulation has been characterized by the prob- 
ability distribution function of the density and pressure 
as well as by computing the velocity and density power- 
spectra and the energy spectrum. However, although very 
informative, these quantities are not easily obtained from 
the observations (and if so, would rely on various assump- 
tions). On the other hand, column density and average 
velocity along the line of sight are more straightforwardly 
inferred from observational signals. Thus, we attempt here 
to further characterize the flow structure by giving some 
statistics of the column density and the velocity along the 
lines of sight. Since our experimental setup is strongly 
anisotropic (the flow enters from the left and right box 
faces and leaves from the top and bottom faces), we al- 
ways show the results obtained along the x-axis and along 
the y-axis. 

2.1. Line of sight characteristics 

Top panel of Fig. Q] shows the column density, N, along 
the y-axis as a function of the x-coordinate (full line) as 
well as the column density along the x-axis as a function 
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Fig. 2. Top panel: average velocity along the line of sight. 
Bottom panel: velocity dispersion along the line of sight. 
Full lines represent lines of sight parallel to y-axis whereas 
dotted lines represent lines of sight parallel to x-axis. 



of the y-coordinate (dotted line). All the lines of sight 
available in the simulation, i.e. 10000, are displayed. The 
average value is about 2 x 10 20 cm -2 . These values are 
comparable to several values quoted in the literature for 
high latitude lines of sight (e.g. Heiles & Troland 2003) and 
indicate that the total amount of gas in the simulations is 
comparable to the total amount of gas along these lines of 
sight. We note, however, that only 20-30 % of the lines of 
sight quoted by Heiles & Troland have such low column 
density. Our simulation is therefore representative of the 
most diffuse lines of sight. 

Bottom panel of Fig. Q] shows the integrated quantity 
/ n/Tdl along x and y-axis. This quantity which traces 
mainly the CNM is directly obtained by HI 21cm line ab- 
sorption spectra. An estimate of the CNM column density 
can be directly obtained by multiplying this quantity by 
the CNM temperature which is about 50 K. Again the 
values obtained are very similar to the values quoted by 
Heiles & Troland (2003) for the lines of sight having total 
HI column densities of about 2 x 10 20 cm~ 2 . 

Top panel of Fig. [2] shows the average velocity along 
the line of sight defined as: 



< v > - 



Ep' 



(i) 



whereas bottom panel of Fig. [2] shows the velocity disper- 
sion along the line of sight: 



< 5v 2 >= 



J2p( v - v o) 2 



(2) 



The mean velocity is close to zero (0.05 km/s for lines of 
sight parallel to y-axis and 0.39 km/s for lines of sight 
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Fig. 3. Histogram of column density. The full lines repre- 
sent the total column density along x and y-axis whereas 
the dotted lines represent the column density (along x and 
y-axis) of CNM only. 

parallel to x-axis) although significant fluctuations (up to 
2 km/s) are present, the standard deviation is 0.9 and 0.7 
km/s for lines of sight parallel to y and x-axis, respectively. 
The velocity dispersion along the lines of sight is signifi- 
cantly larger. The mean value is about 2.29 and 2.5 km/s 
along y and x-axis respectively, whereas the standard vari- 
ation is about 0.7 km/s in both cases. Interestingly, adopt- 
ing a CNM sound speed of 0.7 km/s, these values corre- 
spond to a Mach number of about M. = 2.5 km/s / 0.7 
km/s ~ 3.5 which is similar to the values obtained for the 
CNM clouds by Crovisier (1981) and Heiles & Troland 
(2003). 



2.2. Line of sight statistics 

Here, we present statistics of lines of sight. Figure [3] dis- 
plays the histogram of total column density along x and 
y-axis (full lines) and of CNM (defined as gas having tem- 
perature lower than 200 K) column density (dotted lines). 
Although our simulation described a coherent region of 
20 pc, whereas the lines of sight of the Millennium survey 
are not correlated, it is worth doing some comparison with 
Fig. 7 of Heiles & Troland (2003). Their first panel shows 
that CNM column density peaks at ~ 1.5 x 10 20 cm -2 close 
to our case (~lx 10 20 cm" 2 ). We find that about 20% of 
lines of sight have a small (< 2 x 10 19 cm" 2 ) column den- 
sity of CNM, whereas Heiles & Troland find about 30 % 
of such lines of sight (first panel of their Fig. 7). They find 
that the total HI column density peaks at about 2 x 10 20 
cm -2 , which is close to what we find as well. They find no 
line of sight with column density smaller than ~ 1 x 10 20 
cm~ 2 . In our case the smallest column density is about 
7 x 10 19 cm -2 . Note that the smallest column density of 
atomic gas observed in the Galaxy is about 4.5 x 10 19 
cm~ 2 (Lockman et al. 1986). 

Figure 2] shows the histogram of the fraction, r of 
WNM (full lines), CNM (dotted lines), thermally un- 
stable gas (dashed lines) along the lines of sight, e.g. 
r cnm = n CiCnm /n Ci tot- Note the definition follows those of 
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Fig. 4. Histogram of WNM (full line), CNM (dotted line) 
and thermally unstable gas (dashed line) column density 
fraction. Top panel is for lines of sight parallel to y-axis 
whereas bottom panel is for lines of sight parallel to x-axis. 
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Fig. 5. Fraction of CNM column density as a function of 
total column density. 



paper I, i.e. WNM is gas having temperature higher than 
5000 K whereas unstable gas has temperature between 200 
and 5000 K. Top (Bottom) panel is for line of sight parallel 
to y-axis (x-axis). The fraction of CNM can be compared 
with bottom panel of Fig. 7 from Heiles & Troland (2003). 
Again both appear to be similar. The median value in our 
case is about 0.34, whereas Heiles & Troland quote 0.3. 
We have about 10% of lines of sight for which the fraction 
of CNM is smaller than 0.05 and they find about 30% of 
lines of sight with no CNM. We have about 30% of lines 
of sight with a fraction of CNM larger than 0.5, whereas 
they find this fraction to be about 10%. 

Figure [5] displays the fraction of CNM, r cnm , as a func- 
tion of the total column density along the line of sight. It 
can be compared with Fig. 8 of Heiles & Troland (2003). 
A similar trend is observed in both cases, r cnm increases 
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Fig. 6. Fraction of lines of sight which contain a gas parti- 
cle of a given gas density (upper panel) and mean column 
density along the line of sight (bottom panel) of gas hav- 
ing a given density (only lines of sight which contain a 
fluid particle of the corresponding density are taken into 
account). Full lines represent lines of sight parallel to the 
y-axis whereas dotted lines represent lines of sight parallel 
to x-axis. 



with the total column density. However, the observations 
appear to have a dispersion larger than the simulation re- 
sults. This could possibly be a consequence of the fact that 
the observations probe several regions which are physically 
uncorrelated. 

Altogether, these numbers appear to be similar. One 
should however keep in mind that our statistics are based 
on correlated lines of sight, whereas the lines of sight ob- 
served in the Millennium survey represent regions with a 
priori different physical conditions. As emphasized in pa- 
per I, physical conditions like external triggering and level 
of turbulence do influence the physical parameters of the 
flow as the CNM/WNM fractions. Indeed, since the region 
that we are simulating is actively forming CNM structures, 
it is unsurprising that the number of lines of sight with 
no CNM is underestimated and the fraction of lines of 
sight dominated by CNM is overestimated. Nevertheless, 
the relatively good agreement suggests that the numerical 
setup we used leads to a reasonably realistic medium. 

Another useful statistics is the number of lines of sight 
which cross a piece of fluid of a given density, n. The re- 
sult is shown on top panel of Fig. [5] (logarithmic intervals 
of density are used). While most lines of sight intercept 
fluid particles of density between 0.5 and 100 cm' 3 , there 
is about 10% lines of sight which cross gas denser than 
1000 cm -3 and less than 0.1% which cross gas denser 
than 10 4 cm -3 . In particular, these results suggest that 
the so-called TSAS (Heiles 1997) are somehow rare events 
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as suggested by Stanimirovfc et al. (2003) and Johnstone 
et al. (2003), even if quantitative prediction should be con- 
sidered with great care. 

Finally, an interesting question is: when a line of sight 
crosses a fluid particle of a given density, what is the mean 
column density of gas at this density along the line of 
sight ? The result is shown in the second panel of Fig. [5] 
(note that we have plotted the mean column density of 
gas at given density, po along lines of sight which do 
intercept a fluid particle of density such that logpo and 
logpo + dlogp ). While the mean column densities of gas 
at density 1-2 cm -3 and ~100 cm -3 are approximately 
equal to about 3xl0 19 cm -2 , the mean column densities 
of gas denser than 10 3 cm -3 are about 2 to 3 times larger. 
This means that, although few lines of sight intercept fluid 
particles denser than 10 3 cm' 3 , when it is the case, the 
column density of the dense gas is usually a significant 
fraction of the total column density along the line of sight. 

3. HI spectra 

3.1. Emission and absorption 

Since HI spectra provide most of the available information 
on the interstellar atomic hydrogen, we give here various 
examples of synthetic HI spectra calculated from the nu- 
merical simulation, in emission and in absorption. This is 
certainly the most straightforward way of comparing ob- 
servations and simulations. Figure [7] shows the density and 
the velocity fields along 4 lines of sight parallel to the x- 
axis (only half of the box simulation is displayed since the 
other half contains only WNM and no CNM structure) as 
well as the 21cm line emission and extinction. 

The first case (top panels) shows a line of sight with 
low density CNM structures (~ 30 cm -3 ). Two groups of 
clouds can be seen (x ~ 5.5 and ~ 6 pc), each of them 
being strongly substructured. The relative velocity of the 
two groups of clouds is about 5 km/s and each of them 
presents an internal velocity dispersion of a few km/s. In 
spite of this complexity, the emission spectrum appears to 
be relatively smooth. This is a consequence of the fact that 
the thermal broadening of the HI 21cm line emission is not 
small compared to the turbulent broadening. The large 
component of width ~ 20 km/s is the WNM emission. The 
two narrow components seen in emission and in absorption 
are due to the contribution of the two groups of clouds. 
The largest extinction is about 0.007 and the width of the 
lines is typically around 3-5 km/s. These numbers are very 
similar to the numbers quoted in Braun & Kanekar (2005) 
and Stanimirovic & Heiles (2005). 

The second case (second line of panels) shows a line 
of sight which contains four CNM structures of density 
higher than 100 cm~ 3 . Whereas three of them (x ~ 2.5, 
x ~ 3.2 and x ~ 5.5 pc) present a modest velocity disper- 
sion of few km/s with respect to the others, one of them 
(located at x ~ 5 pc) moves at a velocity of about 4 km/s 
with respect to the other. As a consequence, the corre- 
sponding HI spectra present two peaks, the first one (v ~ 



1 km/s) being due to the contribution of the three struc- 
tures moving at approximately the same velocities and 
the second (v ~ 4 km/s) being due to the fourth structure 
(x ~ 5 pc). Whereas the first peak is slightly broadened 
by the velocity dispersion of the 3 structures, the width 
of the second one is mainly thermal. 

The third case shows about 10 CNM structures denser 
than 100 cm -3 and about the same number of CNM struc- 
tures having a lower density. A velocity gradient of about 
1 km/s/pc takes place between x ~ 1 and x ~ 6 pc. As a 
consequence, the HI spectra present a single peak, instead 
of 2 in the preceding cases, which is significantly broad- 
ened by the velocity dispersion of the structures. Note that 
this spectrum is qualitatively similar (temperature, ex- 
tinction and shape) to the spectra shown by Heiles (2001) 
except for the width which is 2 times larger than the case 
presented here. 

The fourth case shows a line of sight which crosses 
the highest density reached in Fig. 1 of paper II (x ~ 
1.65 pc and y ~ 7.07pc). The HI 21cm line is strongly 
saturated because of the large column density. The spectra 
are again broadened by the velocity dispersion between 
CNM structures. Again, the complex structure of the flow 
cannot be easily recovered from the HI spectra which are 
relatively smooth. 

Figure [8] shows four lines of sight along the y- 
axis. Indeed, since our experimental setup is strongly 
anisotropic, it is worth to see what are the consequences 
of this anisotropy on the synthetic HI spectra. In particu- 
lar, along the y-axis few uncorrelated groups of structures 
can contribute to the synthetic signals. This is well illus- 
trated by the first, second and third cases in which 3 to 
5 groups of structures are present. Each group presents 
an internal velocity dispersion of a few km/s and an av- 
erage velocity with respect to the other groups which can 
be as high as ~ 5 km/s. Interestingly, this is somehow 
similar to the structures obtained by Hennebelle & Passot 
(2006) in the case where significant Alfven waves are ini- 
tially present in the flow. Whereas in their case, this was 
due to the effect of Alfven waves which tend to subfrag- 
mcnt the CNM structures, here the subfragmentation of 
the group of structures is likely due to the turbulent fluc- 
tuations. As a consequence, the synthetic spectra are on 
average broader than the spectra obtained along the x-axis 
(factor ~ 1.3 - 1.5). 

The fourth case (Fig. |5J) corresponds to the line of sight 
parallel to the y-axis which crosses the highest density 
reached in the simulation. As for the fourth case (Fig. [7]), 
the 21cm line is strongly saturated. The narrow peak seen 
in emission is produced by the strong density fluctuation 
(located at x ~ 7 pc) which is also displayed in Fig. 5 of 
paper II. As explained in paper II, we believe that this 
structure is a good candidate to explain the TSAS ob- 
served in the atomic gas. This spectrum could therefore 
constitute a typical signature of such event. Note that in 
the fourth case of Fig. [7j this narrow component is not 
seen because in front of it, stands a less dense structure 
(~ 10 3 cm -3 ) which has a high extinction. On the con- 




Fig. 7. Density and velocity fields along 4 lines of sight (parallel to x-axis) and synthetic HI spectra in emission and 
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Fig. 8. Density and velocity fields along 4 lines of sight (parallel to y-axis) and synthetic HI spectra in emission 
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trary, in the present case, there is no dense structure on 
the right side of the densest peak. 

3.2. Effect of finite resolution 

Here, we investigate the effect of finite resolution on the HI 
emission spectra. We convolve the HI spectra calculated 
along the lines of sight by Gaussian functions of various 
widths at half maximum, b. Figure [9] shows the results. 
Four series of spectra showing the influence of a finite res- 
olution are displayed. The full line of first panel of each 
series shows the same spectra as those shown in Fig. [8j 
(left-top series of spectra is first case of Fig. [FJ left-bottom 
is second case, right-top is third case, whereas left-bottom 
is fourth case). The full line of second, third and fourth 
panels of each series shows the spectra obtained by con- 
volving the spectrum distribution by a Gaussian function 
of width at half maximum equal to b = 10, 20 and 40 re- 
spectively. The dotted (dashed) lines of each panel show 
the unsmoothed spectra shifted from the central position 
by A =-5 (5), -10 (10), -20 (20) and -40 (40) cells respec- 
tively. 

Left-top series (first case of Fig. [3]) show that there 
is a significant variability of the HI 21cm emission spec- 
tra when different lines of sight are selected (dotted and 
dashed lines). This has severe consequences on the con- 
volved spectra which vary significantly when the width 
at half maximum of the Gaussian function used for the 
convolution, increases. As can be seen, the shape of the 
unconvolved spectra (first panel) is rather different from 
the shape obtained with b — 10, 20 and 40. 

Left-bottom series (second case of Fig. [5]) and right-top 
series (third case of Fig. [8j show cases for which the vari- 
ability of the spectra is much less important than in the 
previous case. The general shape of the convolved spectra 
(full lines) does not change much, the fluctuations appear 
to be on average smaller than ~ 10% with some stronger 
fluctuations of about ~ 20 — 30%. The difference between 
this behaviour and what has been found in the previous 
case, is largely due to the fact that the CNM structures are 
much smaller in the latter than in the former. Therefore 
the spatial length over which the column density varies 
significantly, is also smaller in left-top case than in these 
2 cases, making the dependency on the lobe effect more 
important. 

Right-bottom series (fourth case of Fig. [5]) show also 
significant variability in spite of the large emission value. 
This is due to the fact that the line of sight (A — 0) crosses 
a small scale and dense structure (seen in Fig. 5 of paper 
II) which results from a strong collision between 2 CNM 
structures. Therefore, as indicated in Fig. [6j the column 
density varies significantly along the line of sight in few 
cells. 

On average, the smoothed spectra are not much 
broader than the unsmoothed one. This is certainly a con- 
sequence of the fact that the velocity dispersion along the 
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Fig. 10. Mass of the structures as a function of their size. 
One finds the relation M oc L 7 with 7 ~ 1.7. 

lines of sight presents variations which are smaller than 
its average value (Sect. 2.1). 

4. Clouds properties 

In paper II, the mass spectrum of the structures formed 
in the simulation has been measured and theoretical argu- 
ments to explain it, have been presented. Here, we further 
study the structures properties. First, we consider the case 
of structures simply defined by the clipping procedure ex- 
plained in paper I and II. Second in order to quantify the 
fact that the structures are spatially correlated and not 
randomly distributed, we investigate some properties of 
groups of structures. 

4.1. Properties of individual clouds 

We study various structure properties which could be com- 
pared with observations or with other theoretical works. 
We pay special attention to the distribution of column 
density within clouds, since an interesting observational 
result has been obtained by Heiles & Troland (2005) with 
which comparison is possible. 

In the following, we define the cloud size, L, as L = 
y/Ii/M, where M is the cloud mass whereas I\ is the 
highest inertial momentum (I2 being the smallest), i.e. 
the highest eigenvalue of the inertial matrix, X, defined as 
Ixx — J y 2 dm, T yy — J x 2 dm and I xy = T yx = — J xydm. 

4.1.1. Mass-size relation and size distribution 

An interesting property for cloud characterization is the 
mass-size relation. Figure [TU] shows the mass of the struc- 
ture as a function of their size. A clear correlation is seen, 
leading to the relation M oc L 1 with 7 ~ 1.7. 

Since, as will be seen below, this relation turns out 
to be useful to understand the column density distri- 
bution, it is worth to understand its physical origin. 
For this purpose, we consider the distribution function, 
M{L), of structures of size L. The number of structures 
of size between L and L + dL is equal to the num- 
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Fig. 11. Aspect ratio of the clouds as a function of their 
mass. The solid line shows the average value per logarith- 
mic interval of mass. 
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Fig. 12. Aspect ratio of the densest regions as a function 
of the largest gas density within the region. The solid line 
shows the average value per density interval. 

ber of structures of mass between M and M + dM, 
N(L)dL = M{M)dM where M oc LP '. In paper II, we 
obtain that Af (M) oc M~@ with /3 ~ 1.7. Therefore, we 
have: N{L) = N{M) dM/dL oc £-07+7-1. with the val- 
ues of (3 and 7 given above, we have: — /S7 + 7— 1 — —2.19. 
This implies that M(L)L 2 ~ L~ - 19 , i.e. the number of 
structures of a given scale L within a (bidimensional) vol- 
ume L 2 , is almost the same as the number of structures of 
scale L 1 within a volume L' 2 . This implies that the flow is 
nearly scale-invariant (see Elmegreen 1997 and Padoan & 
Nordlund 2002). This is likely a consequence of the flow 
being turbulent, although in the present problem, there 
are characteristic scales which may explain the deviation 
from the scale-invariant behaviour. Another possibility is 
that the size of the simulation is still insufficient and the 
numerical mesh, as well as the injection scale, still have 
an influence. 



4.1.2. Cloud shape 

Figure QT] shows the aspect ratio of the cloud defined as 
\J I2 / Ii ■ The solid line represents the mean value per log- 
arithmic mass interval. The mean aspect ratio is approx- 



Fig. 13. Number of lines of sight (crossing a single CNM 
structure) per logarithmic column density interval. Full 
line corresponds to the (10 4 ) 2 cells simulation, dotted line 
is for the (5000) 2 one, dashed line is for the (2500) 2 cells 
simulation. 

imately constant and equal to about 0.4 until M ~ 0.1 
M s /pc. It then increases to about 0.6 for larger mass. 
Nevertheless we note that there is a considerable disper- 
sion, some structures being very elongated. 

Another related question is the shape of the very dense 
regions observed in the simulations (see section 3.2 of pa- 
per II). Since this high density component is created by 
shocks, it is expected that these regions should be very 
elongated. Figure [12] shows the aspect ratio of these re- 
gions as a function of the peak density (note that since 
these regions are rare, 6 time steps have been used to 
provide this plot) and confirms that these structures are 
elongated with a mean aspect ratio of about ~ 0.3. We 
note that a significant fraction of them has an aspect ra- 
tio smaller than 0.2, which is significantly smaller than the 
mean aspect ratio of unshocked CNM structures. Indeed 
very few CNM clouds have an aspect ratio smaller than 
0.2. This suggests that this could constitute an interest- 
ing test for the assumption of TSAS being shocked CNM 
structures. 



4.2. Column density distribution of clouds 

In the quest for comparison between observations and 
simulations, an interesting result obtained by Heiles & 
Troland (2005), appears to deserve special attention. 
Heiles & Troland show that the column density distribu- 
tion of their CNM structures follows the relation J\f(N) oc 
1/Af. In this section, we study in some details the column 
density distribution of the CNM structures and propose a 
theoretical explanation for this relation. 



4.2.1. Numerical results 

Here, we display the column density distribution, M(N) 
obtained by putting together the column density distribu- 
tion of all the structures identified in our simulation. The 
result is shown in Fig. [13] which shows the number of lines 
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of sight, rii, per logarithmic interval of column density, 
dni /dlogN = Ndni / dN = NAf(N) , for various numerical 
resolutions. 

The low column density part is obviously spoiled by 
numerical resolution effects. In the intermediate part, say 
between 3 x 10 18 cm~ 2 and ~ 3 x 10 19 cm -2 , the distri- 
bution does not appear to depend on the numerical res- 
olution for l g > 5000. This suggests that in this range of 
column density, numerical convergence has been reached. 
On the contrary, the high column density part strongly de- 
pends on the numerical resolution. In particular the higher 
the resolution, the larger the number of lines of sight hav- 
ing high column density. We therefore believe that only 
the intermediate column density part is physically signif- 
icant. 

As can be seen from the highest resolution simulation 
(solid line) displayed in Fig. [I3l the column density distri- 
bution is stiffer between 1.5 x 10 19 and 3 x 10 19 cm~ 2 than 
between 1.5 x 10 18 and 1.5 x 10 19 cm -2 . Since there is no 
clear physical reason for this to be true, this is most likely 
a numerical effect due to the finite size of our numerical 
experiment. Since the size of the simulation box is only 20 
pc and the WNM density about 0.8 cm~ 3 , it is indeed dif- 
ficult to produce structures with column density equal to, 
or higher, than about 20x3.08xl0 18 x0.8 ~ 5xl0 19 cm" 2 . 
Therefore one expects that the number of big structures 
and therefore the number of lines of sight with large col- 
umn density drop above a certain threshold. The value of 
this threshold for the column density is certainly smaller 
than ~ 5 x 10 19 cm~ 2 . This is in good agreement with 
what has been found in paper II for the mass spectrum 
(Figure 15) which presents a cutoff for high mass as well. 

Therefore, we believe that the most reliable part of the 
column density distribution is the one between ~ 1.5x 10 18 
and 1.5 x 10 19 cm -2 . Note that this range of column den- 
sity corresponds to the part of the mass distribution, (Fig. 
15-18 of paper II) for which numerical convergence has 
been reached. The solid line shows that between 1.5 x 10 18 
and 1.5 x 10 19 cm" 2 , one has about Af(N) oc 1/AT 1 ' 15 ~ 
1/7V 12 which is slightly stiffer but nevertheless close to 
the AT(N) oc 1/N distribution. Although the disagreement 
appears to be statistically significant (Heiles 2006, Private 
Communication), it appears to be mainly due to the col- 
umn densities larger than 10 20 cm" 2 , which we cannot 
probe here. It is also possible that the column density dis- 
tribution is biased because of large column density compo- 
nents being blends of two smaller components or because 
of small column density being missed (Heiles 2006, Private 
Communication) . 

Note that Heiles & Troland (2005) infer JV(N) ocl/N 
for N between 2 x 10 18 and 2 x 10 20 cm . The reason of 
this disagreement is probably due, at least in part, to the 
numerical setup and most likely to the finite spatial range 
of our simulation box. Future studies using larger box size 
and larger numerical resolution, may give answer to these 
questions. 
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Fig. 14. Average column density of the structures as a 
function of their mass. One finds the relation < N >oc M v 
with r) ~ 0.52. 

4.2.2. Physical explanation 

Here, we propose a physical explanation for the column 
density distribution, Af(N) oc N~ a . 

For this purpose, we make the assumption that the 
clouds are sufficiently uniform for us to consider that the 
column density within the cloud does not strongly fluctu- 
ate. Figure fT4l displays the mean structure column density 
as a function of their mass, M. A very good agreement is 
obtained with the relation < A" >oc M 5 . 

This result can be understood as follows. Let us con- 
sider a spherical CNM cloud of uniform density n. Its mass 
can be written M ~ nL 2 and the column density N ~ nL. 
Thus < A" >oc M v with r\ = 0.5. In the case we are con- 
sidering, M oc L 7 which implies that the CNM structures 
are not spherical clouds with uniform density. To fix idea, 
we consider 2 cases. 

First, the CNM clouds are roughly roundish but the 
density is not uniform because the internal structure is 
complex. The number of lines of sight crossing the struc- 
ture is proportional to L, thus the mean column density 
follows < N ><x L 1 ~ 1 oc Af^-ri/T. With 7 = 1.7, we 
have (7 — l)/7 ~ 0.41. 

Second, the CNM clouds are filaments with uniform 
density. The major axis has a length L, whereas the mi- 
nor axis has a length L 1 ^ 1 . If the filament is seen along 
the minor axis, the mean column density again follows 
< N >oc L 7_1 oc M (t- 1 )/ 7 . If the filament is seen along 
the major axis, then < A" >oc L oc Af 1 / 7 with I/7 ~ 0.59. 

Since all the values found are slightly above or slightly 
below 0.5, it is not surprising to find rj ~ 0.5. 

The number of lines of sight crossing a structure of 
column density between N and N + dN is proportional 
to the number of structures having a mass between M 
and M + dM times the length of the structures (since the 
number of line of sight crossing a structure of size L is 
proportional to L). Thus we have: Af(N)dN oc Af(M) x 
Af 1 / 7 x dM. This leads to N(N) ~ jV(-0+V'y+i->7)/>7. 
Since with (3=1.7, 7=1.7 and 77=0. 5 we have — (— /3+1/7+ 
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1 - rf)/n ^ 1.22, we get: J\f(N) oc 1/N 1 - 22 . This value is 
remarkably close to what is measured in the simulation. 

The good agreement between the analytical value and 
the exponent measured in the simulation, suggests the va- 
lidity of the approach. However, since we used the value 
7, inferred from the simulation, it could be that its value 
is different in the ISM (as previously discussed, the size of 
the computational box or the numerical resolution could 
be insufficient). In particular, it could be, that real atomic 
flows are closer from strict scale invariance than what 
has been found in the simulation. If instead of taking 
the value of 7 measured in the simulation, we assume 
N(L) oc L~ 2 , thus L~ 2 dL = M^^dM which leads to 
M oc L 7 , 7 = l/(/3 — 1). With this value of 7, we have 
— (— /3 + I/7 + 1 — T))/t] = 1 which is exactly the value in- 
ferred by Heiles & Troland (2005). Interestingly, the result 
appears to be independent of (3 and 77. 

Let us stress that the result, Af(N) oc N~ a , with 
a ~ 1 — 1.2, is largely a consequence of N(L) cx L~ 2 
and of the fact that the column density across a structure 
is nearly proportional to its size, which is due to the ther- 
mally bistable nature of the flow, since in that case, the 
density is determined almost entirely by the mean pres- 
sure and does not depend on scale or mass. Future works, 
using larger simulations, may clarify to which extent such 
flow follows Af(L) cx L~ 2 . 
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Fig. 16. Measure of the internal velocity dispersion of 
structures extracted from numerical data obtained with 
different smoothing and clipping threshold. 



4.2.3. The 3D case 

Since the simulation and the theoretical derivation are 
bidimcnsional, it is worth investigating the value of the 
exponent of the column density distribution in the 3D 
case. The theory developed in paper II for the structure 
mass spectrum, Af(M), predicts that Af(M) oc M~@ with 
(3 ~ 16/9. Assuming scale invariance, we get Af(L) oc 
1/L 3 , which leads to: L~ 3 dL = M~@dM and thus to 
M oc L 7 , with 7 = 2/(j3 — 1). We again make the as- 
sumption that < N >oc M v . For a spheroidal cloud with 
uniform density, we have obviously r/ — 1/3. In 3D, we 
have Af(N)dN = N{M) x L 2 x dM. Thus we obtain 
JV(N) ~ AK-/3+2/7+i-'))Ai = AT 1 , which is again tne 

value inferred by Heiles & Troland (2005) (as in the 2D 
case, it does not depend on (3 or on 77). 

Note that for (3 = 1.8, we get 7 = 2/0.8 = 2.5, implying 
M oc L 2,5 . This is very similar for the mass-size relation 
inferred by Larson (1981) for the molecular clouds. Even 
more interesting is the comparison with the results ob- 
tained by Heithausen et al. (1998) for CO clumps of mass 
ranging from 10 4 solar mass down to Jupiter mass. For 
these clouds, they measure M oc L 23 , Af(M) oc M -1 - 8 
and M[L) oc L~ 3 . This is consistent with what we infer 
for the CNM clouds. 

4.3. Group of structures properties 

The aim of this section is to further characterize the struc- 
ture of the flow and in particular the fact that the CNM 



structures appear to be highly correlated and not ran- 
domly distributed. 

4.3.1. Methods 

In order to achieve this, we smooth the density field by 
convolving it with Gaussian of various widths at half maxi- 
mum, b. This has the advantage of mimicking observations 
done with a telescope having a gaussian lobe and a spa- 
tial resolution equal to b x dx, where dx is the size of 
one of our cells. For b = 40, this leads to a resolution of 
about 0.08 pc. At 100 pc of distance, this corresponds to 
a resolution of 3 arc minutes which is comparable to the 
resolution of the Arecibo telescope. Figure [15] shows the 
density field displayed in Fig. 1 of paper II which has been 
convolved by a Gaussian of width at half maximum equal 
to b = 40 cells. As expected, the field displayed in Fig. [TS] 
is more uniform than the field displayed in Fig. 1 of paper 
II, although significant large scale fluctuations persist. As 
in paper II, we identify the structures by applying a sim- 
ple clipping algorithm, i.e. the structures are defined as 
groups of connex cells denser than a density threshold po. 
Two values of po are considered in the following, namely 
5 and 30 cm -3 . 

We stress that the analysis presented in this section, 
ought to characterize an ensemble of individual, discon- 
nected from each other, CNM structures which are spa- 
tially correlated and would be seen with a telescope of 
insufficient spatial resolution as a coherent cloud. 




Fig. 15. Smoothed density field. The original density field has been convolved with a Gaussian having a width at half 
maximum equal to b — 40 cells leading to a spatial resolution of about 0.1 pc. 



4.3.2. Velocity dispersion 

The cloud velocity dispersion is an interesting parame- 
ter for the understanding and the characterization of the 
turbulence which takes place within the atomic medium. 
Before we pursue, we would like to stress that the use of 
the word turbulence in referring to the gas motions should 
be made with some care. Indeed, the structure of the flow 
studied here is very different from the incompressible flu- 
ids which have been extensively studied and which consti- 
tute our basic conception of the turbulent processes. As 
discussed above, the turbulence is produced by individual 
cloudlets motion rather than by eddies in a nearly homo- 
geneous medium. In particular, the coupling between the 
various spatial scales is presumably strongly modified by 
the stiff density fronts. 

The cloud velocity dispersion is defined as stated by 
Eq. [2l where vq is the mean cloud velocity and where all 



cells having a density above the threshold and lying in 
the area defined by the smoothed density field have been 
counted (we stress that we use unsmoothed data to com- 
pute the velocity dispersion) . Figure [16] shows results for 
b = and 40 and clipping thresholds 5 and 30 cm -3 . 

In the case with no smoothing, we find that most of the 
structures have a velocity dispersion lower than 1 km/s 
meaning that most of structures have a subsonic veloc- 
ity dispersion. By performing a simple linear fit, we infer: 
\/5v 2 ~ 7 km/s x(r/lpc) 8 . Note that the size of struc- 
tures found with no smoothing is generally smaller than 
0.1 pc so that this relation is valid only up to V Sv 2 ^ 1 
km/s. 

With b = 40, the structures are bigger and the veloc- 
ity dispersion higher. As can be seen, it is not strongly 
affected by the clipping threshold (factor ~2). Typical 
velocity dispersions are now around 1 km/s and up to 
3 km/s for the biggest structures. We infer V Sv 2 ~ 1.4 
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Fig. 17. Measure of the cloud filling factor in structures 
extracted from numerical data obtained with different 
smoothing and clipping threshold. 



km/s x (r/lpc) 34 for a clipping threshold of 30 cm 3 
and \/5v 2 ~ 2.2 km/s x (r/lpc) 4 for a clipping thresh- 
old of 5 cm -3 . We note that these values are reminiscent of 
the values observed in HI clouds by Crovisier et al. (1981) 
and Heiles & Troland (2003, 2005), although Figure 12 
of Heiles & Troland (2003) suggests that the mean value 
of the Mach number inside structures is about 3. In our 
case, it is slightly lower, 1.5-2.5, depending on the thresh- 
old which is used to define the structures. At this stage, 
it seems difficult to do a more quantitative comparison 
since the present simulations are only 2D. In particular, 
the approach used in this section, excludes any projection 
effect along the line of sight. 

The fact that the internal velocity dispersion of CNM 
structures is low (< 1 km/s) but that the velocity disper- 
sion of the groups of structures is comparable to observa- 
tions, suggests that the turbulence inside CNM clouds is 
largely due to individual long-living cloudlet motions. 

Again, we note that the velocity dispersion-size rela- 
tion is reminiscent of what is observed in the case of molec- 
ular clouds (Larson 1981). 

4.3.3. Filling factor 

One of the striking aspects of the structures displayed in 
Fig. 1 and 2 of paper II is that the CNM appears to be very 
fragmented even in the densest parts of the flow. In order 
to quantify this effect, we have computed the filling factor 
of the structures extracted from the smoothed data. This 
is obtained by simply computing within these structures, 
the filling factor of the gas (taken from unsmoothed data) 
having a density above the clipping threshold. The result 
is displayed in Fig. [17] which shows the filling factor for 



various b and clipping thresholds. As expected, the value of 
the filling factor presents a significant scatter and depends 
on 6, the width at half maximum of the Gaussian used to 
compute the smoothing. Typical values range nevertheless 
between 0.8 and 0.2 whereas the average value is about 
~ 0.5. 

5. Conclusion 

In a companion paper, high resolution numerical simula- 
tions of the turbulent interstellar atomic hydrogen have 
been presented. Because of the high spatial resolution 
(0.002 pc) reached in these simulations, scales approach- 
ing the size of the smallest structures observed in HI can 
be described. 

In this paper we have further analyzed these numeri- 
cal data in order to make preliminary comparisons with 
HI observations. We provide statistical informations on 
column density, average velocity and velocity fluctuations 
along the line of sight and, when possible, we compare 
with the results of the Millennium Arecibo 21 centime- 
ter absorption line carried out by Heiles & Troland (2003, 
2005). 

We compute HI 21cm line spectra in emission and in 
absorption along eight lines of sight. We reach the conclu- 
sion that although these spectra are relatively smoothed 
the lines of sight from which they are calculated appear 
to be highly complex and fragmented. In particular, the 
"turbulence" is produced by the individual motions of 
few long living CNM structures rather than by motions 
within nearly homogeneous or isothermal medium. This 
suggests that the Mach number deduced from observation 
is more indicative of the velocity dispersion from cloud to 
cloud rather than of the internal velocity dispersion of the 
clouds. We also consider the influence of a Gaussian lobe 
on the HI spectra and find that, it has in some cases, a 
drastic influence. 

To further characterize the CNM structures produced 
in the simulation, we have studied some of their proper- 
ties, namely mass-size relation and shape. In an attempt 
to describe the high level of spatial correlations between 
the CNM structures and to mimic observations done with 
a limited spatial resolution, we have smoothed the data 
simulations by convolving them with a Gaussian func- 
tion of various sizes at half maximum and analyzed the 
results. By doing this, it is possible to identify groups 
of structures that would be seen as a single object and 
study their properties like filling factor and velocity dis- 
persion. We find the mass-size relation, in the bidimen- 
sional simulation, M cx L 1 with 7 ~ 1.7 and based on 
the assumption N{L) cx L~ 3 , we speculate that in 3D, 
we may have M cx L 25 . We also find that the velocity 
dispersion increases with the size of the structures, L, as 
v 7 < Sv 2 > cx L aA . These behaviour are reminiscent of the 
Larson laws's (Larson 1981, Heithausen et al. 1998) in- 
ferred for the molecular clouds. Moreover, since the mass 
spectrum that we obtain for the CNM structures in pa- 
per II, turns out to be compatible with the mass spec- 
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trum observed for the molecular clouds (Kramer et al. 
1998, Heithausen et al. 1998), this suggests that (or at 
least is compatible with it), the origin of the molecular 
clouds, is determined at a very early stage in the diffuse 
atomic medium. This is also consistent with the propo- 
sition made by Falgarone & Puget (1986) and recently 
reinvestigated by Hennebelle & Inutsuka (2006) of molec- 
ular clouds being 2-phase objects. Note that it may be at 
first surprising to propose a theory for the origin of molec- 
ular clouds which does not include self-gravity. However, 
as revealed by CO observations (e.g. Bertoldi & McKee 
1992, Heyer et al. 2001), molecular clouds of mass smaller 
than 10 3 solar mass, are not gravitationally bound mak- 
ing very hard to invoke gravity as an important physical 
agent for the formation of these clouds. Heyer et al. (2001) 
also speculate that an external pressure in needed to con- 
fine them. External pressure is naturally provided by the 
2-phasc physics. 

Finally, we have investigated the distribution of the 
cloud column density, Af(N). We find that between N ~ 
2 x 10 18 cm" 2 and N ~ 2 x 10 19 cm" 2 , Af(N) oc A/"" 1 - 2 
whereas from their survey, Heiles & Troland (2005) in- 
ferred Af(N) oc A^ 1 for N ~ 2 x 10 18 cm~ 2 and N ~ 
2 x 10 20 cm~ 2 . We speculate that the discrepancy in the 
range of column density is due to the setup of the numer- 
ical simulation which limits the highest column densities. 
We also propose that the index of the power law is due 
to the flow being nearly scale-invariant, i.e. N{L) oc L~ D 
(D being the dimension), and the CNM structures being 
roughly uniform in density. Whereas the first is likely to be 
a consequence of turbulence, the second is a consequence 
of the flow being thermally bistable, that is to say the den- 
sity of the CNM does not change with scale or structure 
mass. 
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